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■<^j- ■ 1. Introduction 

l> 

Three years ago, Geim, Novoselov and co-workers have succeeded in isolating and 
contacting a single layer of graphite (graphene) . [I] Contrary to common wisdom, this 
experiment showed that true two-dimensional lattices are thermodynamically stable. 
[21 Ej This stability comes about because the system gently crumples to the third 
direction forming ripples. [I] It is therefore an important problem to study the 
effect of out-plane phonons on the electronic properties of the system. Contrary to 
two-dimensional electron-systems in semiconductor heterostructures, where (scalar) 
electrons interact with bulk[5] or surface[6] phonons, in graphene there are two- 
dimensional in-plane as well as out-of-plane vibrational modes to which Dirac (spinor) 
fermions will couple. 

There are two different vertex-types modeling the interaction of electrons with 
optical phonons. First, due to the atomic displacement within the plane, the tunneling- 
matrix element between two carbon atoms varies. This gives rise to a Su-Schrieffer- 
Heeger type coupling (current-current coupling) . [7] The two-dimensional gauge field is 
composed by the longitudinal optical (LO) and transverse optical (TO) phonons which 
are degenerate at q — 0. This vertex type is also found from symmetry arguments. [8] 



Effect of Holstein phonons on the electronic properties of graphene 



2 



There are also out-of-plane vibrations. These lattice displacements are symmetric with 
respect to their neighboring atoms and thus couple to the electronic densities. The 
optical (ZO) modes can then be described within the Holstein model. [9] For a recent 
account on the Green's function of the Holstein polaron, see Ref. [10j and references 
therein. 

Out-of-plane modes are energetically smaller than in-plane vibrations due to the 
hybridization of the sp 2 -orbitals within the graphene sheet. A tight-binding calculation 
including nearest and second-nearest neighbors yields uzo ~ ^lo/2 for the optical 
branch close to the r-point.[TT] For a first-principle calculation of the phonon-spectra, 
see the work by Wirtz and Rubio.[T2] The effect of the electron-phonon coupling on the 
local density of states of zig-zag graphene ribbons has been studied by Sasaki et al. .[13] 
The effect of the optical phonons on the Raman spectrum of disordered graphene was 
studied by Castro Neto and Guinea. [H] 

The effect of Holstein phonons on transport properties has been receiving renewed 
interest in the context of the one-dimensional Holstein-Hubbard model. [15] Here, we will 
discuss Holstein phonons in graphene for the following reason. It is currently believed 
that transport properties can be well described within a semi-classical Boltzmann 
approach. [T6J [T7J [18] This implies well-defined quasi-particles, i.e., only one band is 
considered and inter-band transitions are ruled out. Doing so, several scattering 
mechanisms have been discussed, ranging from local defects (vacancies or substitutions), 
long-ranged Coulomb impurities in the substrate or due to adsorbed atoms to acoustical 
phonons. P2J [20] 

Optical phonons cannot be treated within the one-band Boltzmann approach [21] 
since they would induce interband-transitions at typical densities of n<5 x 10 12 cm~ 2 . 
It is therefore crucial to assess this scattering mechanism via the Kubo-formalism. In 
section II, we present our model for Holstein phonons and calculate the electronic self- 
energy within the first Born approximation in section III. We further compute the 
optical conductivity in section IV, using the full Greens function, but neglecting vertex 
corrections. We close with conclusions and outlook. 

2. The effective model 

In two-dimensional graphene sheets, first-principle calculations reveal that the Born- 
Oppenheimer approximation is not valid for doped graphene sheets. [22] Frohlich 
polarons[23j are thus not a good starting point to describe electron-lattice interaction 
in graphene. Here, we present the study of the electron-lattice coupling due to localized 
Holstein-phonons in a two-dimensional honeycomb lattice, thus treating the ZO-phonons 
(out-of-plane vibrations) . 

The honeycomb lattice (the lattice of graphene) is made of two interpenetrating 
triangular lattices, defining two non-equivalent sites, usually labeled as A and B sites. [23] 
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The model Hamiltonian for ZO-polarons in graphene reads as follows: 

H-tY, [al{Ri)b a {Ri + 8)+ H.c] + Y"A C * + V e P (!) 

i,crS ' Q 

where a) a (Ri) (6j.(J2j)) creates an electron at an atom of the A(B) sub- lattice and c q 
are creation phonon operators. The energy u q is the dispersion of the ZO-phonon and 
V ep the electron-phonon interaction. We model the coupling to the ZO-phonons as the 
usually density-density coupling [25] 

V ep = DY, [4(fliK0RO + tf a {Ri)b a (Ri) - l] Qj , (2) 

with Qj defined as 

Qj = YX q e iqBd (c q + cl q ), (3) 

and 

(4) 



\ 2MNu q 

and M the ion's mass and iV the number of unit cells in the crystal. Note that we 
couple the density with respect to the half-filled band such that particle-hole symmetry 
is conserved by reversing the sign of the coupling constant. In the following, though, we 
shall neglect this shift in the bosonic operators. 

The effect of Holstein polarons is obtained using oj q ~ c^o and transforming the 
operators in Hamiltonian (JT]) to momentum space. This gives 

H = - t Y. (<f><Pq,* b W + H ' C -) + E ^>C q C q + Vep , (5) 
q,a q 

with (p q = J2s e~ tS q , 8 the vectors connecting the three nearest neighbors on the 
honeycomb lattice, and V ep given by 

V ep = D X q (a^ap+q^ + b p a b p+q ^ (c q + cl q ) . (6) 

Let us comment on the coupling constant D. Due to the mirror symmetry of the 
graphene-sheet, one might think that the linear coupling to lattice displacements is zero. 
But the mirror symmetry is broken for samples where graphene lies on top of a SiC>2- 
or SiC-substrate (for a discussion, see Ref. [8]). To quantify the coupling constant 
in terms of the dimensionless constant g = y/NDXo/ujQ, we assume that the coupling 
mechanism is due to a variation of the hopping matrix element. This yields g of the 
order of unity. [1] 



3. Second order perturbation theory 

If the phonon energy scale is much smaller than the electronic energy scale, Migdal's 
theorem states that it is sufficient to calculate the lowest-order self-energy diagram. [26] 
This diagram can further be calculated using the bare electron Greens function. 
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Still, the importance of electron-phonon coupling also depends on the dimensionality 
of the system. E.g., in self-assembled quantum dots, vertex "corrections" to the 
polarization of an electron- hole pair give rise to charge-cancellation. [27] thus changing 
the optical conductivity significantly. [28] Also for A 3 C G0 (A=K,Rb)[29j and for general 
one-dimensional systems. |30j Migdal's theorem is not valid. Furthermore, electron- 
electron interaction can affect the effective electron-phonon interaction. |15j 

In graphene sheets, electron-electron interaction is generally neglected, i.e., one 
assumes a "normal" ground-state at zero doping (one electron per unit cell) - 
characterized by a semi-metal. Since the average kinetic and interaction energy per 
particle both scale with s/n where n is the carrier density, the interaction does 
not become important at finite doping, either. Electron-electron interaction is also 
neglected in recent works on localization [31] even though disorder enhances the effect 
of interaction. [32] The same should hold for the electron-phonon vertex corrections 
which results in an effective electron-electron interaction and we thus believe that the 
assumption of Migdal's theory is a good starting point to discuss the Holstein-phonons 
in graphene. 

Because of the existence of two sub-lattices, the Green's function needs to be written 
as a 2 x 2 matrix: 

r (h \-( G AA,a( k i T ) G AB ^(k,T) \ 

a[ ,T) \G BA Ak,r) G BB>a (k,r) ) ' [() 

with 

GaaA^t) = - (Ta fciCT (r)4, CT (0)> , 
G AB Ak,r)= - (Ta fc)(T (r)6^(0)) , 
G BA Ak,r) = - (T& fc , CT (r)4 j(T (0)> , 

G BB Ak,r)= - (Tb hta (T)b{ >a {0)) , (8) 
where r is the "imaginary" time, and T is the time ordering operator. 

Up to second order in perturbation theory, and after transforming the time 
dependence of the Green's function to Matsubara frequencies, we obtain the following 
result: 

G = G(°' + GW, (9) 
where the matrix S is defined as 

s = ( l AA l AB ) , (10) 

\ ^BA ^BB j 

and the matrix element S Q/ g is defined by (a, (3 = A,B) 

S^(^ n ,p) = -^D 2 X^(°)(q, W )GS(p-q,^-^). (11) 

Because both X q and D^(q,w) are momentum independent, i.e., we set uj q = uq 
and the phonon propagator D^(q, iv) is given by 

° m <«» = 7-5^73, (12) 
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the matrix elements of the self energy matrix can be written in a simplified form, reading 
T, a p(iu> n ) = —jjJ2 q ,u D 2 'XqD^ (iiy)G^(q,iuj n +iu) where X is X q with u q replaced 
by ojq. Notice that we neglected the constant shift in the self-energy due to the shift of 
the bosonic displacement operator, still present in Eq. (TIT]) . This is consistent since we 
also neglected the Hartree correction to the self energy. 

From the above we can write down a Dyson equation for the electronic propagator, 
given by 

G = G (0) + G (0) SG, (13) 

which has to be solved for G. The equation giving the Matsubara Green's function for 
the free electronic system reads 

/ iu n -t<p{k) 
\ —t(f)*(k) iuj n 

and 1 is the 2x2 unit matrix. Within the Dirac cone approximation, which applies 
in a range of 1 eV, one has T, AB (iu n ) = T, BA (iu n ) = 0. Also the following result holds 
T, AA (iu n ) = ^BB^n) = E(iu; n ) leading to a simplified form for the Dyson equation 
ffl3|) . which can be readily solved, giving 

/ iu n - S(zw n ) t<p(k) \ 
G _ V t(p*(k) iu n - E(^ n ) ) ^ 
(iu n - E)(iu n - £) - t 2 \<p\ 2 



G®{iu n ,k) = l, (14) 



The matrix elements of Eq. (Tl5l) can be put into a simpler form reading 

1/2 

G AA (u n ,k) = " — x ..i xi , (16) 

jt±l lU} n ~ M^n) -Jt\(j){k)\ 

G AB (u n ,k) = £ -■ °. J .+ U(hM , (17) 

jt^l «*>n - -jt\(j)(k)\ 
j e -i5(k) j2 

G BA {uj n ,k) = Y ' "°' 



iu n - E(iuj n ) -jt\(f>(k)\ ' 
G BB (uj n , k) = G AA (u n , k). (19) 

The summation over the bosonic frequency v in Eq. can be performed using 
standard methods leading to 

,jr±i 2 U*. -#10(4)1 +"0 iw„-j*W<l)|-'V7 

The integrals over the momentum variable in Eq. (|20|) can be easily computed at zero 
temperature and within the Dirac cone approximation, yielding an explicit form for the 
self energy. 



Effect of Holstein phonons on the electronic properties of graphene 



6 



3.1. Zero doping 



At zero temperature and zero chemical potential (that is at the neutrality point) the 
self energy, denoted by £ — > So, has a simplified form reading 

1 1 



, Q [lUJ ri 



Nh 2 y 2 



lUJr, 



t 



+ 



r 



[q)\/h - Uq/% iuj n + t\(j)(q)\/h + uj 

where we have introduced the missing h's omitted in the beginning of this section and 
g denotes a dimensionless coupling constant of order unity, defined below Eq. (JBJ). 
Performing the analytical continuation iuj n — > uj + i5 and computing the momentum 
integral in Eq. ( 1211 one obtains the retarded self energy, E ret (o;), of the polaron problem, 
which reads 



iret , 



UJ 



Ac (9^o\ 2 (_^_ 



+ 



2tt V h ) 

UJ + UJq/% 



vi- 



lli 



(v F k c 




VK- 



2tt V h 



(22) 



6{uj — uj /h)9(vFk c + ujq/Ti — to) 
uj — uj /h)9(vFk c + ujq/Ti + uj) 



a 2 3-y/3/2, and a 



1.42 A. 



with Tik 2 c = (2tt) 2 /4, A c 

One can easily obtain the renormalization of the electronic spectrum due to the 
retarded self-energy induced by the phonons, using Rayleigh-Schrodinger perturbation 
theory. [25] In this scheme, the uj dependence of the self-energy is replace by the bare 
electronic dispersion. Close to the Dirac point, the dispersion is simply given by 
uj = ±vpk where the upper (lower) sign holds for electrons (holes) and we obtain the 
new energy spectrum E(k) from 

E{k) = v F hk + h^E T et {v F k) (23) 

for the conduction band and 



E(k) = -v F hk + M?E r ct (-v F fc) 



(24) 



for the valence band. Electron-hole symmetry is thus preserved, if 3?Eq 



ret i 



-3?E r ct 



x) which is the case (see Eq. (|22|) ). 
Considering that vpk is the smallest energy in the problem, that is k is very close 

re 




to the Dirac point, we obtain for 9ftEQ Ct (t>^fc) the result 



TXVi 



2 , v F k c 
In 



h J uJa/fi 



(25) 



Considering the out-of-plane optical mode of graphene [12J, which has an energy of 
ujzo — 0.1 eV and using g 2 = 10 we obtain 



A, 



TTVf 



h 



2 v F k c 
In ■ 



uj /h 



0.02 . 



(26) 




Figure 1. (Color online) Renormalized energy band at the neutrality point, given by 
E(k) = ztvphk ± Ti^tH^^Fk) for two different values of the coupling constant g. 



In Fig. dj the renormalized energy dispersion is shown for various values of the coupling 
constant. We also included a large coupling constant g = 10 for better illustration of 
the effect of the electron-phonon interaction. 

3.2. Finite doping 

At finite doping, the one-particle dispersion j\<fi(k)\ must be replaced by j\(f>(k)\ — \x in 
Eqs. ( TT6i) - ( |T9l) and Eq. ( |20l) . where n = TivpkF denotes the Fermi energy. The retarded 



self-energy can then be written as £ 
have 

AW A Ac f 9UJo 



ret i 



UJ 



yret 
^0 



u /h 



In 



uj + /i/fr) + AE(w + fi/h) where we 



2vr V h 

(n/h — uj — uo/h)(u — ujo/h) 



(uj /ny 



a, 



%TX- 



2tt V h 



(fl/% - UJ + uj /h)(uj + UJq/K) 

2 / !,) I,l„ \ 

9(uj + uj /h)9(fi/h 



UJ 



UJQ_ 




ojq/Ti — uj) 



uj — uj /h)9(fi/h + uj /h — uj) 



(27) 



We note that for finite doping (uj) diverges logarithmically at u — ±uq/H. 



4. Optical conductivity 

In this section we want to compute the optical conductivity of graphene and study 
how the conductivity of the system is affected by the out-of-plane (ZO) phonons. To 
determine the conductivity one needs to know the current operator j x , which is composed 
of the paramagnetic and diamagnetic contributions j x = j£ + A x (t)j^, each of them 
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3* = E - 3)4(fc)M*0 " (<f(*0 - 3)6i(fc)o.(fc)] , (28) 

fc,(T 

and 

jf = £ + 3)4(fc)6.(fc) + (0*(fc) + 3)6+ (*K(*)] . (29) 

The Kubo formula for the conductivity is given by[34j 

xx[ } iA s {u + 10+) + ihA s {u + «0+) ' 1 j 

with A s = N C A C the area of the sample, and A c the area of the unit cell, from which it 
follows that 

9f£<x(^) = Dd(w) + ^— ^ , (31) 

where D is the charge stiffness which reads 

D = ~ 71 — n hX s ' (3) 

The incoherent contribution to the conductivity A xx (u + i0 + ) is obtained from A xx (iuj n ), 
with this latter quantity defined as 

= f dre^ < T t3 p x (r)jf(O) > . (33) 



o 



The relevant quantity 5A ra (w + i0 + ) is given by 

%a xx (u + to + ) = t -^4- (4)] E 

[n F {e) - n F (e + uj%)} A Xl (k, u + e/h)A X2 {k, e/h)f(k, A 1; A 2 ) , (34) 
with A x (k,uj) given by 

A x (k,cu) = -2SG x R (k,u + iO + ) , (35) 
G x (k,iu) given by 

G A (fc, z^ n ) = — — r — ) r-pr , (36) 

ioj n - £(za; n J - (Xt\(p{k) \ - n)/n 

and f(k, A 1; A 2 ) given by 



/(fc, A 1? A 2 ) = 2|^(fc) - 3| 2 - AiA 2 
Using the fact that 



(fc) - 3)2 ^ + ( ^ ) - 3)2 W T) 



E>(fcMW*oi) = Y,<p*(k)g(\m\) = 1 -Y,\m\ 2 g(\m\) & 

k k 6 k 
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Figure 2. (Color online) Optical conductivity for two different values of the coupling 
constant g at zero doping (T = T = 10~ 4 eV). Inset: The relative conductivity 
Act = a(g — 10) — a(g — 0) with and without $RE(a;). 



where is an arbitrary function of the absolute value of </>(&), and the fact that 

in the Dirac cone approximation one has 

2 (fc) 



J2tt/3 



Ik) 



cos(2(9) -isin(20)] 



(39) 



and a similar expression for the complex conjugate expression [0*(fc)] 2 /|0(fc)| 2 , one 
obtains a simplified expression for (l34"j) . given by 



t 2 e 2 a 2 f de 



16ti 2 



2nh 



E E 



fc Ai,A 2 =±l 

[n F (e) - n F (e + u)%)} A Xl (k, u + e/h)A x ' 2 (k, e/h) 
[18-2|0(fc)| 2 (A 1 A 2 + l)]. 



(40) 



The expression ( |4*0|) is valid only in the Dirac cone approximation, and therefore 
one must replace t\(p(k)\ by vphk and the integral over the momentum can be easily 
done. If one ignores the effect of phonons the calculation is straightforward, leading to 
a conductivity of the form (at zero doping) [33] 

" 2 < (hw) 2 \ , ( ojh ^ 



Kct(w) 



7r e 



1 



tanh 



(41) 



2 /i ^ ist 2 y V 4fc s T , 

One should note that from Eq. (T4TI) one has 9far(0) = for finite T. This is not seen in 
Fig. [2] because the temperature scale is very small. 

The solution of Eq. PU|) allows us to determine the optical conductivity taking into 
account the effect of Holstein phonons. In addition, we mimic the effect of impurities 
by adding a small imaginary part T to the self-energy 

9E(w) = 3£ ep M - T , (42) 

where 3S ep (u;) is obtained from Eq. ( J22l) . 
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There are two types of momentum integrals in (1401) . which have the form 

rkc k 2n+1 




(43) 



with n = 0,1. The analytical expressions are given in the appendix. The final energy 
integration can be done numerically. 

In Fig. El the conductivity is shown for various coupling constants g at zero doping 
and low temperature T and damping T due to impurity scattering. There is a drop 
of the conductivity (relative to the conductivity of a clean system) starting at around 
the phonon energy u and reaching a constant value for u>>2u} . Without the real part 
of the electronic self-energy, there is a pronounced peak at twice the phonon energy. 
This is shown in the inset of Fig. El where we plot Act = a(g = 10) — o~(g = 0) 
for T = T = O.OOOleV and including, respectively, not including 3££(a>) in the above 
expressions. It is thus crucial to include the full self energy in the renormalization of 
the particle Green's function. 

We thus obtain as main result that there remains no pronounced structure in 
the conductivity due to ZO-phonons if the full self-energy is used. We attribute this 
fact to the apparently asymmetric way the self-energy enters in the Green's function 
with respect to the electron (j=l) and hole (j=-l) channel which destroys possible 
interferences between the two carriers (see Eq. (1221) ). 

In Figs. [3]and|H the optical conductivity is shown for different values of the coupling 
constant g at \i — 0.05eV and \i = O.leV, respectively. The insets show the relative 
conductivity due to the electron-phonon interaction, Act = a(g = 10) — a(g = 0), 
including, respectively, not including 9ft£(u;) in the above expressions. Again we see a 
distinct difference in the result due to the renormalization of the electron-hole spectra. 

It is clear that the results at finite chemical potential are markedly different from the 
results at the neutrality point. At finite chemical potential the system is characterized by 
a Drude like behavior followed by a strong increase of the conductivity when the photon 
frequency reach the value of twice the chemical potential. At zero doping there is no 
Drude weight and the system response is characterized only by inter-band transitions. 
We see only weak renormalization of the Drude peak due to the Holstein phonons as 
well as negligible effects at finite frequencies. We note that the results for dto~(u) when 
g = were first obtained by Peres et al. [24j and by Gusynin et aL|35j 

5. Summary 

In this paper, we have calculated the effect of Holstein polarons on the electronic 
properties of graphene. Holstein polarons arise through the coupling of out-of- 
plane optical (ZO) modes to the conduction electrons, described as Dirac Fermions. 
Throughout this work, we assumed Migdal's theorem to be valid and calculated the 
electronic self-energy within the first Born approximation. We find that the Fermi 
velocity becomes renormalized within 1 %. 
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Figure 3. (Color online) Optical conductivity for two different values of the coupling 
constant g at fj, = 0.05eV (T = T = 10~ 4 cV. Inset: The relative conductivity 
Act = a(g = 10) — a(g = 0) with and without 5RS(w). 
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Figure 4. (Color online) Optical conductivity for two different values of the coupling 
constant g at fj, = O.leV (T = T = 10 _4 eV. Inset: The relative conductivity 
Act = a(g = 10) — a(g = 0) with and without 5RS(a;). 



The main purpose of this work was to assess the effect of Holstein phonons on 
the conductivity within the Kubo-formula. Due to the large phonon-energy, electron 
scattering from Holstein phonons induces interband transitions for usual carrier densities 
(corresponding to a gate voltage of ~ 50V) and can thus not be treating within the 
one-band Boltzmann approach. We thus calculated the optical conductivity within the 
Kubo-formula, employing the full Green's function but neglecting vertex corrections. We 
find a pronounced kink-like peak at twice the ZO-phonon energy if only the imaginary 
part of the self-energy is considered. This peak vanishes when the real part of the self- 
energy is included. Further, we see only weak renormalization of the Drude peak due 
to the Holstein phonons as well as negligible effects at finite frequencies. We conclude 
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that scattering from Holstein-phonons can be neglected in usual transport properties. 

The effect of lattice vibrations on the electronic properties of graphene is still 
not fully understood. Especially the coupling of substrate phonons [36] or in-plane 
oscillations to the conduction electrons is interesting. In the later case, this will lead to 
a non-diagonal electronic self-energy due to the current-current coupling. 
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Appendix A. Momentum integrals 



The momentum integrals have the form (n = 0, 1) 

k 2n+1 dk 



[(A — k) 2 + B 2 } [(C — k) 2 + D 2 ] 
The general solution yields 



(A.1) 



In ~ G 



F? tan" 1 



A - fc 



+ F 9 n tan" 1 



C-fc* 



- F 3 n In ((A - k) 2 + B 2 ) - F'l In ((C - k) 2 + D 

where the denominator reads G = (B 2 + (A - C) 2 ) 2 + 2 ((A - C) 2 
the factors are given by 

AD{A 3 - 2A 2 C - 2B 2 C + A[B 2 + C 2 + D 2 )) , 
AB((B 2 + (A - C) 2 )C + (-2A + C)D 2 ) , 



(A.2) 

B 2 )D 2 + D 4 and 



Fl 



F 4 ° = 2BD(-A 2 -B 2 + C 2 



D 



and 



Fl 
F\ 
F\ 
F\ 



AD(A 5 - 2A A C + 2B A C + A 3 (2B 2 + C 2 + D 2 ) + AB 2 (B 2 - 3(C 2 + D 2 ))) , 
AB((B 2 + (A - C) 2 )C 3 + C(-3(A 2 + B 2 ) + 2C 2 )D 2 + (2 A + C)D 4 ) , 
2BD(A 4 - AA 3 C - AAB 2 C + B 2 (B 2 - C 2 - D 2 ) + A 2 {2B 2 + 3{C 2 + D 2 ))) 
2BD{C 2 (3(A 2 + B 2 ) - AAC + C 2 ) - {A 2 + B 2 + AAC - 2C 2 )D 2 + D 4 ) . 
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